arXiv:1506.01795v3 [gr-qc] 6 Mar 2016 


Higher Spins Tunneling from a Time Dependent and 
Spherically Symmetric Black Hole 


Haryanto M. SiahaarEl 

Physics Department, Parahyangan Catholic University, 

Jalan Ciumbuleuit 94, Bandung 40141, Indonesia 

Abstract 

The discussions of Hawking radiation via tunneling method have been performed ex¬ 
tensively in the case of scalar particles. Moreover, there are also several works in discussing 
the tunneling method for Hawking radiation by using higher spins, e.g. neutrino, photon, 
and gravitino, in the background of static black holes. Interestingly, it is found that the 
Hawking temperature for static black holes using the higher spins particles has no differ¬ 
ence compared to the one computed using scalars. In this paper, we study the Hawking 
radiation for a spherically symmetric and time dependent black holes using the tunneling 
of Dirac particles, photon, and gravitino. We find that the obtained Hawking temperature 
is similar to the one derived in the tunneling method by using scalars. 
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1 Introduction 


General relativity predicts the existence of black holes, i.e. objects from which nothing can 
escape. Hawking found that incorporating quantum mechanics changes this understanding 
about black holes mm- He showed that black holes are not completely black. They radiate 
thermal spectrum due to the quantum effects. Furthermore, Hawking found a connection 
between the surface gravity of a black hole and its corresponding temperature. 

There are several ways in understanding the mechanism of black holes to radiate. One 
of the latest approaches was given by Parikh and Wilczek[31 Sj, where the radiation of black 
holes is described as a quantum tunneling effect of particles through the event horizon. This 
method is found to be simpler and more intuitive. In the Parikh-Wilczek (PW) method, which 
is called sometime as the radial null geodesic method, one first computes the across horizon 
tunneling amplitude as an exponentiation of the imaginary part of the corresponding particle’s 
action in the outgoing mode. Then the principle of detailed balance is used to connect the 
tunneling amplitude with the Boltzmann factor, hence the temperature can be obtained. 

In fact, the radial null geodesic method is not the only way to describe the Hawking 
radiation by using the tunneling mechanism. There is an alternative, known as the Hamilton- 
Jacobi method [5], where one solves the semiclassical equation of motion for the tunneled 
particles by using the Hamilton-Jacobi ansatz for the corresponding particle’s wave function. 
The discussions for Hamilton-Jacobi method in discussing the Hawking temperature and black 
hole thermodynamics for various cases and black holes have been performed in [g B m El da 
m- However, most of them are confined into the tunneling of scalar fields where one starts 
with a Klein-Gordon equation in curved spacetime and solves it by using Hamilton-Jacobi 
ansatz. 

For static spacetimes, the tunneling mechanism for Hawking radiation has been extended 
to the case of spin ^ fermion [12], as well as photon and gravitino m- These works motivate 
us to study the Hawking mechanism for time dependent black hole’s radiation as the tunneling 
processes of massless higher spins particles. The works presented in this paper can be also 
considered as an extension of our previous work m where the authors study the scalar tun¬ 
neling which gives rise to the Hawking temperature of Vaidya black holes. Since we keep the 
spherically symmetric and time dependent spacetime metric in this paper to be quite general, 
i.e. spacetime whose line element can be written as eq. (12.511 . hence the results presented in 
this paper should be relevant to this class of spacetimcl- 

The organization of paper is as follows. In the next section, we review the null geodesic 
method for a general time dependent metric background. In section jj] we discuss the tunneling 
of massless Dirac particles across the time dependent black hole’s horizon. By using the solu¬ 
tions of massless Dirac fields and then the detailed balance principle, we can get the Hawking 
temperature of the black holes under consideration. Interestingly, the Hawking temperature 
derived in this section is invariant compared to the one obtained in the case of tunneled scalar 
fields |14j . The same prescription is repeated in section B where the starting point is the 
source free Maxwell equation in curved space. Again, the Hamilton-Jacobi ansatz is used to 
solve the corresponding equation. After using the Lorentz gauge condition, we can get the 
solution for the vector fields. The same finding is appeared in section B where the obtained 
Hawking temperature has no difference to that derived in the scalar tunneling M- 

In section [5l the tunneled object is massless gravitino, i.e. particle with spin |. We start 
from the Rarita-Schwinger equation in curved space, and Hamilton-Jacobi ansatz helps us to 

^For example to those discussed in | 16 | . 
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get the solutions. As we would expect after performing the analysis in sections [3] and HI the 
Hawking radiation due to the massless gravitino tunneling yields the same temperature as one 
in the scalar case M- The conclusions are given in the last section. 


2 Radial Null Geodesic Method 


In [3] , Parikh and Wilczek presented a direct and short derivation of Hawking radiation as a 
tunneling process. To dispel the coordinate singularity that the Schwarzschild coordinate has, 
they employ the Painleve transformation, which yields the Schwarzschild spacetime transforms 
to 

ds'^ = — dt^ + 2\j^^^dtdr + , (2.1) 

where dH| = dd^ + sin^ OdcjP' is the metric of 2-sphere with unit radius. In getting the last line 
element from the Schwarzschild metric, the Schwarzschild time tg is transformed as [3] 

= . ( 2 . 2 ) 

For a general spherically symmetric and static spacetime 

ds^ = -A (r) dt^ + Y (r)“^ dr^ + r^d^l , (2.3) 

the Painleve transformation (|2.2I) can be written as 

dt^dt- . (2.4) 

In this paper we work out the tunneling prescription to explain the Hawking radiation for 
a time dependent and spherically symmetric black holes. In general, such black hole solutions 
can be read as HZj 

ds^ = —X {t, r) dt^ + Y {t, r)~^ dr'^ + Z {t, r) dH^ . (2-5) 


In this paper we study a class of time dependent spacetime which has coordinate singularity 
at r which yields Y{t, r) = 0 and X(t, r) = 0. In order to get rid of this coordinate singularity, 
we employ the “generalized” Painleve transformation (|2.4p which is in a differential form can 
be read as 


, dtp dtp 
dtp = -^dt + -^dr . 


( 2 . 6 ) 


For a general spacetime metric (|2.5|1 . it is clear that both dtp/dt and dtp/dr would be some 
functions of t and r. However, by employing the relation which are applied in static case, i.e. 


dtp = dt 


l-Y 

XY 


dr, we find that the mapping 


dtp ^ 
dt 

and 

dtp I 1 — Y {t, r) 
dr y A {t, r) Y {t, r) 


(2.7) 
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( 2 . 8 ) 













remove the coordinate singularity in the general metric (12.51) . As in the static case, the left hand 
side of eq. (12.81) is singular for the vanishing X(t,r) or Y(t,r). This is a normal consequence 
in such transformation to remove the coordinate singularity. 

Furthermore, to get an integrable tp, from the fact that = 0 in eq. (|2.7I1 one realizes 
that the right hand side of eq. (|2.8I) must be independent of time t, sa}!! 

l-Y{t,r)=X{t,r)Y{t,r)C{r) , (2.9) 

where C{r) is an arbitrary function of r. Hence, a set of Painleve transformation (12.7p and 
(12.8p works only for a class of time dependent spherically symmetric spacetime whose metric 
functions satisfy the conditior@ (|2.9|) . 

Consequently, after employing the transformation (12.7p and (12.8p , the metric (12.5p now can 
be written as 


+ + Z {t, r) . 


( 2 . 10 ) 


It is understood that the coordinate singularity at X(t,r) = Y{t,r) = 0 has already been re¬ 
moved in (12.lOp . where instead to be singular, now the gtr component of the line element above 
is indeterminate!. Accordingly, the null radial geodesics from the “Painleve transformed” met¬ 
ric (j2.inp can be written as 


dr 

dt 



( 2 . 11 ) 


and -!-(—) signs denote the outgoing(ingoing) geodesics. Moreover, for a practical benefit one 
can Taylor expand the coefficient X and Y at the near horizon, i.e. 


X{t,r)\^ ~ X'{t,rh)\^ {r - rt) + O (^{r - 


( 2 . 12 ) 


and 


F(t,r)|^ - y'{t,rh)\^ (r - ru) + O (^(r - ^ , 

where is the radius of event horizon. By using the Taylor expansions 
outgoing null radial geodesic (|2.11l) can be approached as 


(2.13) 

(12.1211 and (12.1311 . the 


fl’Y' 1 _ 

2 irh,t)Y' {rh,t) (r - r/,). (2.14) 

Now we use the prescription by Parikh and Wilczek in getting the Hawking temperature 
using the picture of tunneling particle through the event horizon, which is sometime called the 
Radial Null Geodesic or Parikh-Wilczek (PW) method. However, it was found by Chowdhury 

®In the case of static spacetime, eq. (ITOl) is automatically satisfied. 

'*In [Appendix A| we show a constraint for Vaidya black hole mass function m{t,r) which comes from this 
condition. 

®We consider the case where after employing the L’Hospital’s rule, one can get a non singular form of gtr in 
the metric (HIS]). Otherwise, the method developed in this paper might not work since dr/dt in (12.111) could be 
singular. 
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in m that the expression ImS" = f Prdr which appears in original PW method is not canon¬ 
ically invariant. Therefore, we use the prescription by Akhmedov et al [20] in computing the 
tunneling rate which reads 


r ~ exp 


Prdr 

h 


(2.15) 


The term inside of square bracket above then can be computed as 


Im 


Prdr = Im 


Pr H 

j dp'rdr = ‘^dr , 


(2.16) 


where we have made use of the Hamilton equation dr/dt = dH/dpr\r related to the canonical 
variables r and pr (in this case, the radial component of the radius and the momentum). 
Different from the discussions of several authors for a static black hole mass, e.g. Refs. |3] 
and [18j . the outgoing particle’s energy must be time dependent for black holes with varying 
mass. So, the dH' integration at (I2.16|) is for all values of outgoing particle’s energy, say from 
zero to +E (t). 

By using the approximation (12.140 . we can perform the integration (|2.16p . For dr inte¬ 
gration, we can perform a contour integration for the upper half complex plane to avoid the 
coordinate singularity r/j. The result is 


Im 



47r.E (t) 

y/X' {rh,t)Y' {rh,t) 


(2.17) 


Equalizing the tunneling rate (I2.15p with the Boltzmann factor exp [—(dE (t)] for a system with 
time dependence of energy we obtain 


Th 


hy/X' {rh,t) Y' {rh,t) 
Att 


(2.18) 


This temperature is also derived by Nielsen and Yeom in m by using a slightly different 
approach of PW method for a general time dependent background. 

At this point we understand that the temperature (12.181) doesn’t care about the spin of 
tunneled particles. As long as it follows the radial null geodesic, which we know the path 
of massless particles, then it will give the same contribution the temperature measured by 
a detector at infinity independent of the spin that it has. However, the spins that we are 
mentioning here are 0,1/2,1 and 3/2 only, since in the author’s knowledge there is no works 
for Hawking radiation in the tunneling picture which uses spins > 2 as the tunneled particles. 
In the next sections, we will reproduce the temperature (|2.18l) by considering the tunneling of 
Dirac fermion, photon, and gravitino from a spherically symmetric and time dependent black 
hole by using Hamilton-Jacobi method. 


3 Massless Dirac Particle Tunneling 

In this section we will study the Hawking radiation of time dependent black holes where the 
tunneled particle has spin We start by writing an action describing massless Dirac fields in 
curved spacetime [22] 

s^ = I d^xV^mr 'k. (3.1) 
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The corresponding equation of motion can be written as 


= 0 


(3.2) 


where 








(3.3) 


and = I [ 70 , 7 /?]. We use diag(—,+,+,+) as the Minkowski metric tensor and the flat 
spacetime Dirac matrices 7 “ are 


7° = 


72 = 


i 0 
0 -i 


0 

1 0 


a 




7^ = 


0 cj3 

3 0 


<7 


0 Ct2 
2 0 


a 


(3.4) 


The flat Dirac matrices 7 “ and the “curved” ones 7 ^ are related by 7 ^ = 607 “. 

For the general time dependent metric (12.5j) . the tetrads e“ can be expressed as 

\ 


e^ = 


/ Vx 0 00 

0 i/W 0 0 

0 0 r 0 

\ 0 0 0 rsin 0 / 


(3.5) 


where = e'^^e\,riab- Clearly Ca is just the inverse of (|3.5p . Therefore, the Dirac matrices 7 ^ 
using the tetrads (|3.5D can be written as 


f = 


Vx 


/ 1 0 0 
010 
00-1 
V 0 0 0 


0 \ 

0 
0 

-1 / 


(3.6) 


= VF 


7 = 


T = 


r sm ( 


/ 0 0 1 0 \ 

0 0 0 -1 

10 0 0 

V 0 -1 0 0 / 

/ 0 0 0 1 \ 

0 0 10 
0 10 0 

V 1 0 0 0 / 


/ 0 0 0 -1 \ 

0 0 10 

0-100 
V 1 0 0 0 / 


(3.7) 


(3.8) 


(3.9) 


Now we employ the Hamilton-Jacobi ansatz for the spinor wave function describing the 
massless Dirac particles. To simplify the computation, we perform the tunneling process 
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analysis for spin up particle and spin down one separately!^. Explicitly, the Hainilton- 
Jacobi ansatz for our spinors can be read as 




/^\ 

0 

B 

V 0 y 


exp 


h 


S, 


= 


( 0 \ 
c 
0 

\v ) 

where the function Su and are expanded in terms of h as 


exp ( 


Sk — Sok + ^Sik + h S2k + ^ s^k + • • • 


(3.10) 


(3.11) 


(3.12) 


and k is understood as the spin of Dirac particles under consideration, i.e. k = u or k = d. The 
coefficients A, B, C, and V in general are t, r, 0 and (j) dependent, as well as the corresponding 
actions for spins up and down Su and Sd respectively. As we understand in the Dirac formalism, 
(I3.10p is the wave function for spin up fermion, and (13.111) is for the spin down. In the next step, 
we work out the analysis for spin up only, since an analogous method can also be performed 
for spin down which produces the same Hawking temperature. The case of radial null geodesic 
yields the corresponding particle action does not vary with respect to the 6 and cj) coordinates. 
Therefore, by inserting the spin up wave function (|3.ini) into the equation (13.21) . we have a set 
of equations 

^9i5o„ + VFH5,5o„ + 0(fi) = 0, (3.13) 


^fYAdrS^u 


iB 


dtSou + O (h) — 0, 


(3.14) 


We focus only on the leading order terms in the equations, hence 0{h) can be neglected. 
Moreover, the equations (j3.13l) and (I3.14p can be shown in a matrix operation as following 


f iX-^/^dtSou Y-y^drSou = 

V Y-Y^drSou -iX-Y^dtSou ) [ B ) B ) 


(3.15) 


The vanishing of last equation is guaranteed if the determinant of D is 


zero, which leads to 




(3.16) 


We notice that equation (I3.16p is just the equation for scalar particle’s action in curved space 
after we employ the Hamilton-Jacobi ansatz and consider only the leading terms in the equation 
UM- Moreover, the last equation can be rewritten as 

drSou = ±-j^dtSou (3.17) 

where (—)-|- signs correspond to the (outgoing)ingoing modes. The discussions of these modes 
can be found in Apendix 2.A of [23]. 

®Tlie reader sliould be familiar with the Dirac spinor as a direct sum = 'iu © 4'd. 
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In [H], the authors have derived the solution for an equation like (I3.17P where they consider 
the tunneling of scalar particles from a time dependent black hole. Therefore, the techniques 
presented in m can be adopted to get an expression for Su- We look for a general form of 
solution for the actiorfl 

t 

Sou {t,r) = - J E (t') dt' + Sou {t, r) , (3.18) 

0 

where E{t') stands for the time dependent energy of the Dirac particle which tunnels across 
the event horizon. The time dependence of energy is understood since the mass of black hole 
decreases as time passes. Taking the derivative with respect to time in both sides of the last 
equation provides us 

dtSou {t, r) = -E {t) + dtSou {t, r) , 
and from the differentiation with respect to radius r we have 

drSou {t, r) = drSou (t, r) . 

The chain rule allows us to write 

dSpu {t, r) _ dSpu {t,r) ^ dSpu (T r) dt 
dr dr dt dr 

In this section we don’t use the Painleve transformation as we have used in the previous 
section. Therefore, the corresponding radial null geodesic in the background (12.511 is 

^ = ±VJ^. (3.22) 

The +(—) signs in the left hand side of (|3.22l) refers to the geodesic of outgoing(ingoing) null 
particles respectively. Combining equations (I3.2ip and (I3.22p gives us 


(3.19) 

(3.20) 

(3.21) 


dSpu jt, r) _ dSpu {t, r) 1 dSou{t,r) 
dr dr \fXY dt 


(3.23) 


with — (+) signs refer to the outgoing(ingoing) spin up Dirac from the black hole. 

Recall that for the dynamics of a particle with Hamiltonian H and action S, one can show 
the relation Ol] 

dS 

— + H = 0. (3.24) 

The last equation also emerges in the semiclassical discussion, for example, in WKB approxi¬ 
mation to solve the one dimensional Schroedinger equation 



2m dx^ 


+ H(x)T, 


(3.25) 


where we use the ansatz T = and S is the classical action of particle associated to the 
wave function T. From eq. (|3.25p . we may observe that the partial derivative of action with 

discussion of Schrodinger equation with time dependent Hamiltonian which supports this general form 
of action is given in [Appendix B[ 
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respect to time would be a negative quantity for a particle with positive energy, since the 
eigenvalue of H must be positive. Therefore, the (—) sign in (|3.17p belongs to the ontgoing 
particle, 

dtSou = -V^drSou , (3.26) 

since the momentum pr = drSo is positive. Correspondingly, the one with (+) sign refers to 
the ingoing particle, 

dtSou = VJ^drSou . (3.27) 

Then we use equations (|3.2Up . (|3.2ip . (|3.26p . and (|3.27p to get 


dSou {t,r) _ ^E{t) 


dr 


Vjce 


whose solution can be read as 


Sou {t,r) = ±E{t) 


dr 




(3.28) 


(3.29) 


The + and — signs in the last equation belong to the outgoing and ingoing particle respectively. 
Accordingly, a solution for the action (|3.18p can be read as 


Sou 


{t,r) = - [ (t') di'± 


(3.30) 


Plugging the solution (13.3011 into (13.101) gives us 


^uAn — 


and 




u,out — 


f A\ 
0 
B 

V 0 / 

/ A\ 
0 
B 

\ 0 y 




exp 


i 

h 




(3.31) 


(3.32) 


2 

Making the ingoing probability Pin = \'^u,in\ is unity, i.e. all fields that come close to a black 
hole will be absorbed, yields 


J E {t') dt' 


iirE (t) 

\fxEp>' 


Therefore the ontgoing probability can be written as 


Pout — I Obit I — GXp [ 

The “detailed balance” principle tells us that 

Pout = e~^^Pir, 


A-kE (t) \ 

hVxY'J 


(3.33) 


(3.34) 
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which then allow us to write the Hawking temperature for a general time dependent black hole 
(j2.5p as 


Th = 


h^X' {t,rh) y it,rh) 
Att 


(3.35) 


The Hawking temperature (13.351) is interpreted as the measured temperature by a detector at 
infinity where the radiation consists of massless quantum particle with spin ^ moving outward 
radially from the black holes. 


4 Photon Tunneling 

In [13], Majhi and Samanta discuss the tunneling of photon and gravitino which yield the 
Hawking radiation from a static black hole. One of the conclusions in their work is that the 
Hawking radiation in the form of the tunneling of photon and gravitino yields the Hawking 
temperature which is invariant compared to the one computed in the case of scalar tunneling. 
In this section and the next one, we show that the same conclusion is obtained for a time 
dependent black hole. We start from an action for Maxwell fields in curved spacetime, 

S = -\j (4.1) 

Taking the variation of in the action above, we obtain 

V^F^" = 0, (4.2) 

which is known as the Maxwell equation in the absence of the source J'^. By following Majhi 
et al m, we use the Hamilton-Jacobi ansatz for the vector field 

~ ^ (4 3) 

where is the polarization vector. This polarization vector is independent of the spacetime 
coordinates. As usual, the action is expanded as 

OO 

S {t,r,9,(j)) = (4.4) 

i=0 

just like what we did in the spinor case (13.121) . In [13j . the polarization vector is also 
expanded in H since the authors discuss the quantum correction which comes from the higher 
order terms in h of S{t,r,9,4>) and k^. However, since we are not interested in pursuing such 
quantum correction, the polarization vector k^ in (14.31) can be kept np to its semiclassical 
value only. 

Plugging the ansatz (14.3p for the gauge fields into the equation (14.21) which alternatively 
can be expressed as 

d^Ff^u ^ = 0 , (4.5) 

one can get 

{k^d^^So - k>^d''So) df,So = 0. (4.6) 

In getting the last equation, we have taken the limit —)• 0 in the equation (14.51) . We choose 
to work in the Lorentz gauge, 

= 0, (4.7) 
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which after plugging the gauge fields (14.3p gives 


= (4.8) 

Again we have employed the limit ^ 0 in obtaining the equation (14.8|) . In this Lorentz gauge 

condition, the reading of equation (|4.6K reduces to 


fe" (9^5o)(a^5o) = 0. 


(4.9) 


Working on the t and r sectors only in the spacetime under consideration allow us to write 
(14.91) as 

g^\dtSof + g^^idrSof = 0. (4.10) 

We find that the last equation is similar to (j3.17p if we replace Sq with Squ- It is clear since 
the blackground of spacetime where the vector probes come and fall into a black hole is also 
the same, i.e. the metric (|2.5I) . Therefore, the solution for S{t,r,9,4>) for the vector fields 
can be read as 

t 

So {t,r) = - J E {t') dt' ± ■ (4-11) 

0 

Accordingly, the ingoing and outgoing solutions for the vector fields can be read as 


~ exp I - I - / .E {t') dt' — in 


E 




and 


^out ~ 


exp ^ ^ — I E (t') dt' + in 


j 

E{t) \ 


^/xY' j 


(4.12) 


(4.13) 


respectively. The unit incoming probability Pin = the relation between incoming and 

outgoing probabilities 

4nE {t) \ 


Pout — Pin 6 Xp I 




(4.14) 


and the “detailed balance” principle Pout = Pin exp (—/?P(t)) yields the reading of Hawking 
temperature is 


Th = 


h^/X' {t,rh)Y' {t,rh) 


An 


(4.15) 


One observes that the temperature in the last equation is equal to the one computed in the 
Dirac particle case (|3.35p and scalar case 


5 Gravitino tunneling 

We start with the action of massless Rarita-Schwinger fields in curved spacetime, 


/' 


= I d^xVY^Y'' ( df. + 


A>r 


(5.1) 
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where and the Dirac matrices 7 ^ in the action above are those used in the Dirac action 
m- Accordingly, the action (|5.ip tells us that the corresponding equation of motion for 'I'q 
can be read as 

= (5.2) 

which is known as the massless Rarita-Schwinger equation in curved space. It looks like the 
Dirac equation, with the Dirac spinor \k is replaced by the vector-spinor The number of 
degree of freedom of is sixteen, which eight of them are removed by the two additional 
constraints: = 0 and = 0 . 

The Hamilton-jacobi ansatz for vector-spinor can be read as 


and 






( 

0 

\ 0 


( ^ 
C, 
0 


\ 

/ 

\ 

/ 




(5.3) 


(5.4) 


where are the Rarita-Schwinger fields with spins -1-3/2 and —3/2 respectively. 

In the background (j2.5p . equation (j5.2ll for radial geodesic can be read as 


^dtSoiu) + VYB^drSoiu) +Oih)=0, (5.5) 

iB 

VYA^drSoiu) - ^dtSo(u) +Oih) = 0. (5.6) 

The action 5'o(„) is understood as the zeroth order term in the action expansion = 

CXD _ _ 

Moreover, the last two equations are very close to (j3.13D and (|3.14p since the 

i=0 

operator that applies to the massless Rarita-Schwinger fields in (j5.2D is just the same with 
that applies to the massless Dirac field in (j3.2p . Analogous to the technique applies to the 
Dirac fermion in the previous section, we rewrite the equations (|5.5p and (|5.6I) in the form 


V T-V29,5o(„) -iX-V29,5o(,) )\B^ )- \ ) 

whose solution for the action 5o(„) finally can be found as 


Bo{u) (^) ^) 



iirE (t) 

Vx^' 


(5.7) 


Procedure that is needed in obtaining the solution (|5.7p is obvious since we deal the same 
equations in section [3l In discussing gravitino, we replace the complex valued functions A and 


II 













B with the vectors and which are also complex valued. Moreover, the solutions for spin 
I fields then can be written as 




fi{Uyin) 


and 




fl(u,OUt) 


f A, \ 

0 

V 0 y 

/ A, \ 

0 

V 0 y 


exp 


- J E {t') dt' - 


ittE (t) 
y/X’Y' 


exp 


- E {t') dt' + 


iivE (t) 


(5.8) 


(5.9) 


Again, using the ’’detailed balance” principle for the relation of outgoing and ingoing 
probabilities one can get 


Th 


h^/X' {t,rh)Y' {t,rh) 

dvr 


gravitino 


(5.10) 


as the Hawking temperature from dynamical black holes (12.511 in the form of massless gravitino 
tunneling from black holes. We observe that the Hawking temperature due to the massless 
gravitino tunneling is equal to the temperatures computed in the last two sections. 


6 Conclusion and Discussions 

We have analyzed the Hawking radiation in the form of Dirac fermion, photon, and gravitino 
tunneling across the event horizon of time dependent black holes. The resulting Hawking 
temperatures are invariant compared to the one obtained in the scalar tunneling case. The 
results lead to a conclusion that the Hawking temperature obtained in the tunneling method is 
independent of the spins of the tunneled particles. We confirm that for a time dependent and 
spherically symmetric black hole whose metric has the form (12.5p . the Hawking temperature 
is independent of the spins of the tunneled particle, as in the case of static one. 

It is interesting to note that the PW method presented in section [2] works for a limited case 
only, i.e. the spacetime metric whose metric functions satisfy the condition (12.9p . The obtained 
Hawking temperature in this method is confirmed by the result derived via Hamilton-Jacobi 
method by using higher spins tunneling in the succeeding sections as well as scalar tunneling 
[14j . In the Hamilton-Jacobi method we do not make use of the Painleve transformation, hence 
there is no integrability condition (j2.9p to be fulfilled. Such condition also does not appear 
in and yet the same result for Hawking temperature (I2.18P is achieved by the authors. 
Hence, presumably there are some other transformations which yield the line element (12.5p to 
be regular at event horizon, and therefore the PW method can be performed. Clearly each 
of these transformations will have a set of integrability conditions which justify a class of 
spacetime which fits in the computation of Hawking temperature using PW method. Finding 
an alternative to the generalized Painleve transformation and (12.8K which yields the metric 
(j2.5p transforms to be regular at the event horizon would be challenging and we address this 
issue in our future work. 

In their seminal paper [3], Parikh and Wilczek included the back reaction effect in their 
analysis of Hawking radiation in the tunneling picture. Hence, discussing a correction to the 
entropy of time dependent black holes which comes from the back reaction effect should also 
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be possible. For the case of Vaidya black holes entropy, Zhang et al [29] had carried out an 
analysis which takes back reaction effect into account, where the Vaidya spacetime is written 
in the Eddington-Finkelstein coordinate rather that Schwarzschild-like one. We note that 
Zhang et al consider the particle energy to be time independent, unlike the time dependent 
case presented in this paper. We address the analysis of time dependent black hole’s entropy 
by considering back reaction effect in Schwarzschild-like coordinate in our future project. 
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Appendix A Example: Vaidya black holes 


The Vaidya black hole is an interesting object to be studied m- It is an example of exact 
time dependent black hole solution where it solves the Einstein equations with T^i, = pk^k^- 
Since k^k^ = 0, one can consider that the Vaidya solution describes a non empty spacetime 
outside of a null radiating mass. Scalar tunneling as a mechanism for Hawking radiation from 
a Vaidya black hole is also discussed in [26] , and the authors of m study the spin ^ tunneling 
as a Hawking process for this black hole. 

Related to the generalized Painleve transformation (12.71) and (j2.8p . we can take Vaidya 
black holes [251E], 




/m {t, r) 


2 



2m {t, r) \ 
r J 


dt^ 


2m{t,r) 


-1 

-(- r‘^clil.2 ■ 


(A.l) 


as a concrete example. In the last equation, f{t,r) is an arbitrary function that depends on 
the model of mass m{t, r) being used and satisfies 





2m {t, r) \ 
r ) 


dm (t, r) 
dr 


(A.2) 


It follows that from eq. (ESI), the mass functions of Vaidya black hole which are compatible 
with the time dependent Painleve transformation and (1^ should obey an equation 

which reads 


2 m(t,r)/(t,r)^ 


rm {t, r) 


2 


A _ 2m {t,r) 
\ r 


C{r) . 


(A.3) 


Hence, the case of Vaidya black holes whose Hawking temperature can be computed using the 
method proposed in this paper are those with the mass functions that satisfy eq. (IA.3p . which 
is found to be 


Th 


h (m — rm') 
2TTr‘^ 



2 





d {m/ f) r [r — 2m) 


dr 


m — rm' 


1/2 


r=rh{t) 


(A.4) 


Note that at the Schwarzschild limit [25], where /(t, r) = —rh, the Hawking temperature above 
reduces to the familiar Th = for Schwarzschild black hole. 
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Appendix B Time dependent Hamiltonian in Schrodinger equa¬ 
tion 

In this appendix we elaborate some subtleties in time dependent Schrodinger equation in 
supporting the general form of solution for action (I3.18p . The analysis presented here can be 
found in [2H] chapter 31. 

Let us consider a Hamiltonian H[r{t)) which describes a physical system. The variable 
r{t) are parameters which slowly vary with respect to time t. Hence, we understand that 
the Hamiltonian under consideration here is time dependent. Accordingly, the Schrodinger 
equation to be solved in this case is 

. (B.l) 

The associated Hamiltonian operator and energy eigenvalues relation can be read as 


H (r {t)) |/c, r (t)) = Ek (r (t))\k, r (t)) (B.2) 

which is valid at an instantaneous time t only. The state \k, r (t)) is normalized accordingly to 


(/c,r (t) I l,r{t)) = 5ki 

The solution to equation (jB.ip can be expressed as 

^ -i \ dt'Ek(r{t')) 

(t)) =^o.k{t)e 0 \k,r{t)) . 

k 

Plugging (IB.4P into the Schrodinger equation (IB.Ill gives us 

f Q 

k 


a-k + ak-^ \k,r{t)) = 0. 


Using the adiabatic approximation, the solution to Ofc(t) can be read aq 

Uk it) = , 

where Gk (t) is some purely real time dependent function, i.e. 


t 

Gk (t) = i j dt' {k,r \k,r[t')) . 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


(B.7) 


Finally, the solution to Schrodinger equation (jB.ip in the adiabatic approximation can be read 
as 


, , ,, i! dt'{-En{r(t'))+h(k,r{t')\-^\k,r(t'))} 

|T {t)) = e 0 \k,r{t)) , 

which justifies the general form of action’s solution (13.1811 . 


(B.8) 


®A detail derivation can be found in 
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